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I.  INTRODUCTION 


Current  kinetic  energy  (KE)  penetrator  projectile  such  as  the  M829  use  fins  made  of 
aluminum  to  provide  flight  stability.  Typically,  these  fins  are  very  thin  (3.5  mm  is  typical) 
and  consequently  do  not  have  sufficient  mass  to  absorb  and  conduct  away  the  high  heat 
loads  due  to  inbore  and  free  flight  aerodynamic  heating.  The  fins  can  ablate  and  cause 
erratic  flight  due  to  distortion  of  the  fin  or  lack  of  sufficient  fin  area  to  maintain  stability 
and  spin  control. 

The  U-S.  Army  Ballistic  Research  Laboratory  (BRL)  quasi-three-dimensional  heat 
transfer  code1  computationally  models  both  the  inbore  and  free  flight  aerodynamic  heat¬ 
ing.  The  code  provides  designers  of  high  velocity  projectiles  the  capability  to  predict  the 
unsteady  thermal  response  of  the  fins  to  the  flight  environment.  The  code  has  the  ability 
(although  not  yet  perfected)  to  simulate  melting  of  the  fin  with  a  moving  boundary  that 
recedes  as  the  fin  material  reaches  the  melting  temperature. 

Graphs  depicting  the  temperatures  on  the  leading  edge  of  the  fin  and  temperature 
contours  demonstrating  the  heat  transferred  throughout  the  entire  fin  illustrate  the  capa¬ 
bility  of  the  code  and  provide  a  convenient  comparison  for  evaluation  of  the  computational 
results. 

The  boundary  conditions  used  for  this  study  are  described.  Sample  computations  for 
heat  transfer  are  shown  to  illustrate  the  response  of  the  aluminum  fins  to  inbore  and  free 
flight  conditions.  A  comparison  was  made  between  an  aluminum  oxide  coated  fin  and  a 
silicone/fiber  coated  fin.  The  results  demonstrate  that  the  silicon/fiber  coating  provides 
substantial  thermal  protection  compared  to  that  provided  by  the  standard  aluminum  oxide 
coating.  The  results  also  demonstrate  that  the  computational  code  provides  KE  projectile 
designers  a  useful  tool  for  the  evaluation  of  the  unsteady  thermal  response  of  fins  for 
proposed  design  changes  of  fin  and  coating  material.  The  computations  were  carried  out 
using  the  Cray  2  supercomputer  located  at  the  U.S.  Army  Ballistic  Research  Laboratory. 
Aberdeen  Proving  Ground.  MD. 


II.  THEORY 


KE  projectile  fins  are  predominately  made  of  high  strength  aluminum  which  has  the 
physical  properties  of  high  thermal  conductivity  and  relatively  low  melting  point  (950 
K).  These  fins  have  been  computationally  modeled  as  having  either  an  aluminum  oxide 
or  silicon-fiber  coating  which  is  placed  on  the  fin  to  provide  thermal  protection  from  the 
effects  of  inbore  and  aerodynamic  heating. 


1.  INBORE 

The  quasi-three-dimensional  code  was  first  modified  to  simulate  inbore  heating.  A 
flame  temperature  of  3500  I\  was  entered  into  the  code.  This  represents  the  flame  temper¬ 
ature  of  the  propellant  used  in  the  120mm  tank  gun. 
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The  equation  for  determining  the  temperature  at  the  surface  of  the  aluminum  was 
derived  by  combining  the  equations  for  convection  heat  transfer  from  the  boundary  layer 
to  the  coating  and  the  conduction  heat  transfer  from  the  coating  to  the  aluminum. 

T0  =  flame  temperature  of  propellant 

T+  —  temperature  of  coating 

T;  =  temperature  of  aluminum  fin  outer  surface 


convection  conduction 

(t,  -  Ti)  +  (r,+  -  X.-)  =  T0  -  t; 

Q  =  hAAT  convection  heat  transfer 
Q  =  kA^  conduction  heat  transfer 


where: 

Q  =  heat  flux 

A  =  local  area  over  which  heat  transfer  takes  place 
Aj  —  thickness  of  coating 

k  =  thermal  conductivity 

h  =  coefficient  of  convective  heat  transfer. 


Substituting  for  AT 


let 


Q  Q  Ax 

hA  +  kA  ~  °  ‘ 


9.(L  +  *£)=T  _T- 

A\h  +  k  )  *  ‘ 


_  1  Ax 

R=h  +  ~r 


The  quantities  Q,  A,  T0,  Ax,  k,  and  h  are  known  since  the  code  solves  for  Q  using 
information  at  the  last  time  step.  Thus  the  boundary  condition  T“  is  now  established. 

R  is  the  resistance  to  heat  transfer  and  was  computed  using  the  given  coating  thick¬ 
ness,  thermal  conductivity  and  coefficient  of  convective  heat  transfer.  The  heat  transfer 
coefficient  was  determined  from  a  previous  study.2  The  value  used  was: 


2 


h 


G153.6  — 
mi 


J  oulc 
—  6tc  —  K 


The  code  simulated  the  inbore  heating  for  10  milliseconds.  This  is  approximately  the 
amount  of  time  that  the  MS29  projectile  experiences  the  flame  temperature  of  3500  K.  The 
temperature  data  at  10  milliseconds  was  saved  in  a  data  file  and  used  as  initial  conditions 
for  the  free  flight  case. 


2.  FREE  FLIGHT 

The  same  equations  for  determining  the  temperature  at  the  aluminum  surface  were 
used  in  the  free  flight  case.  However.  Tc  was  no  longer  the  flame  temperature  and  the 
coefficient  of  convective  heat  transfer  changed.  T0  became  the  adiabatic  wall  temperature. 
This  value  was  obtained  by  considering  the  results  of  a  three  dimensional  Navier-Stokes 
solution3  for  a  swept  fin  geometry  at  Mach  3.5.  This  solution  also  demonstrated  that 
there  was  very  little  change  in  the  adiabatic  wall  temperature  in  the  chordwise  direction. 
Therefore  a  constant  value  was  assumed.  The  results  were  extrapolated  to  the  velocities  of 
interest  as  determined  from  Firing  Table  FT  120-D-l  (August  19SG).  The  launch  velocity 
at  standard  conditions  is  1G70  m/s. 

The  adiabatic  wall  temperature  changed  with  time  due  to  a  decrease  in  velocity  of 
the  projectile.  Due  to  the  small  time  difference  (1.9  seconds  to  reach  3000  meters),  the 
variation  of  the  adiabatic  wall  temperature  with  time  was  approximated  as  the  following 
linear  function 


T0  =  Tau.  (t)  = 


T  - 

x  au', 


320 

—  x  fane 


where:  Tau>  =  the  value  at  the  gun  muzzle  velocity. 

The  initial  adiabatic  wall  temperature  was  entered  into  this  equation  as  1S30  K.  The 
value  had  decreased  to  1510  I\  at  a  distance  of  3000  meters. 


The  heat  transfer  coefficient  also  changed  with  time  using  the  following  linear  rela¬ 
tionship 

700 

——.x  time 

2.0 


hitag  —  hstagi 


where: 


h,tag  =  convective  heat  transfer  coefficient  at  the  stagnation 
point  on  the  fin  leading  edge  at  an  arbitrary  time  in 
the  trajectory 

h,tagi  =  initial  value  of  the  convective  heat  transfer  coefficient 
at  the  fin  stagnation  point  at  gun  launch,  baaed  on 


3 


previous  studies.2"4  [6000.0  J/(m2  —  sec  —  A')] 
time  =  time  since  gun  launch,  seconds. 

The  coefficient  of  convective  heat  transfer  for  the  trailing  edge  remained  fairly  constant 
and  was  assumed  to  be  constant  based  on  previous  studies.2-4 

The  change  in  the  coefficient  of  convective  heat  transfer  as  a  function  of  chord-wise 
position  was  modeled  using  a  cosine  function  as: 

hloc  —  hmtn  (h max  ^  min  )  ^  ^  (COS  ( Ac  X  7T )  -j-  1.0) 

where: 

hioc  =  local  value  of  convective  heat  transfer  coefficient 
hmxTt  =  minimum  value  of  h  [1800.0  J /{m2  —  sec  —  I\)] 
hmax  =  maximum  value  of  h  (  =  ha(ag) 

Ac  =  length  of  chordwise  distance  for  change  of  h  to  occur,  normalized. 
The  physical  properties  used  for  the  materials  in  this  study  are  given  in  Table  1. 


Table  1.  Physical  Properties. 


Property 

Units 

Al 

AI2  O3 

Sihcon/Fiber 

k 

W/m-K 

173 

12.114 

G.34S9 

P 

kg/m3 

2800 

Cv 

J/kg-K 

869 

Lh 

J/kg 

4518S 

Ax 

m 

0.000127 

0.000254 

where: 

k  =  thermal  conductivity 

p  =  density 

cv  =  specific  heat 

Lh  =  latent  heat  of  melt 

Ax  =  coating  thickness 
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III.  RESULTS 


The  initial  computational  grid  is  shown  in  Figure  1.  The  fin  planform  and  thickness 
are  similar,  but  not  identical,  to  that  of  the  MS29.  The  fins  remained  intact  during 
inbore  heating  with  both  coatings.  The  5  and  10  millisecond  temperature  contour  plots, 
shown  in  Figures  2-5.  clearly  demonstrate  that  much  less  heat  is  transfered  to  the  fin  with 
the  silicon/fiber  coating  than  with  the  aluminum  oxide  coating.  The  leading  edge  of  the 
aluminum  oxide  coated  fin  reached  a  temperature  of  4C0  K  while  the  leading  edge  of  the 
silicon/fiber  coated  fin  only  reached  a  temperature  of  between  310  K  and  340  K. 

The  maximum  effective  range  listed  for  the  M829  is  3000  meters  and  is  reached  in 
approximately  1.9  seconds.  The  free  flight  temperature  contour  plots,  however,  extend 
to  3.0  seconds  to  better  demonstrate  the  effectiveness  of  the  silicon/fiber  coating.  The 
maximum  temperature  of  the  aluminum  fin  was  the  melt  temperature  of  aluminum  which 
was  input  at  950  I\.  The  free  flight  temperature  contour  plots  are  shown  in  Figures  G-17. 
During  free  flight,  the  aluminum  oxide  coated  fin  had  already  melted  a  great  deal  at  1.5 
seconds.  The  silicon/fiber  coated  fin  remained  intact  through  two  seconds.  Very  little,  if 
any.  indication  of  fin  melting  is  apparent  either  at  2.5  seconds  or  3.0  seconds.  It  is  not 
apparent  that  the  fin  had  started  to  melt  at  these  last  two  time  periods  since  the  9th 
contour  had  not  appeared. 

The  temperature  as  a  function  of  time  in  free  flight  at  the  leading  edge  of  the  fin 
located  near  the  midpoint  between  the  fin  root  and  the  fin  tip  is  shown  in  Figure  IS.  The 
leading  edge  of  the  aluminum  oxide  coated  fin  reached  a  temperature  of  460  K  while  the 
leading  edge  of  the  silicon/fiber  coated  fin  only  reached  a  temperature  of  between  310  K 
and  340  Iv  during  the  in-bore  phase  at  the  selected  position.  In  free  flight,  the  silicon/fiber 
coated  fin  reached  a  maximum  temperature  of  620  K  while  the  aluminum  oxide  coated  fin 
reached  a  maximum  temperature  of  S70  K  at  the  selected  position. 


IV.  CONCLUSIONS 


The  results  of  the  numerical  computations  support  the  following  conclusions. 

1.  Melting  did  not  occur  for  the  aluminum  oxide  coated  fin  or  the  silicon/fiber  coated 
fin  while  inbore. 

2.  The  silicon/fiber  coated  fin  did  not  reach  melt  temperature  through  the  weapon 
system’s  maximum  effective  range  while  the  aluminum  oxide  coated  fin  reached  melt  tem¬ 
perature  between  0.5  and  1.0  seconds  of  flight. 

3.  The  silicon/fiber  coating  provides  substantially  enhanced  thermal  protection  com¬ 
pared  to  the  protection  provided  by  the  aluminum  oxide  coating. 
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V.  CLOSING  REMARKS 


The  problem  of  KE  projectile  fin  ablation  is  one  that  will  continue  to  exist  with 
the  improvement  and  redesigning  of  KE  projectiles.  Battlefield  conditions  will  require 
higher  velocity  and  greater  effective  range  of  armor  piercing  rounds.  The  computer  code 
used  in  this  study  can  provide  designers  of  these  new  projectiles  the  ability  to  predict  the 
unsteady  thermal  response  of  stabilizing  fins  to  the  combined  effects  of  inbore  and  free  flight 
environments.  The  code  is  sufficiently  flexible  to  provide  the  capability  to  model  different 
fin  geometries,  construction  materials  and  coating  materials.  The  computational  results 
indicate  that  the  siiicon/fiber  coating  provides  substantially  enhanced  thermal  protection 
compared  to  the  standard  aluminum  oxide  coating.  Although  the  thermal  protection 
provided  by  the  standard  aluminum  oxide  coating  is  shown  to  be  marginal  for  the  flight 
conditions  considered  in  this  computational  study,  it  should  be  pointed  out  that  experience 
has  shown  that  the  aluminum  oxide  coating  does  provide  protection  for  the  aluminum  fins. 
The  protective  qualities  of  the  aluminum  oxide  coating  are  highly  effective  for  protecting 
the  fragile  aluminum  fin  during  the  inbore  phase  of  flight.  The  type  of  protection  provided 
is  believed  to  be  primarily  abrasion  resistence  since  the  fin  is  initially  submerged  in  the 
propellant  bed.  Abrasion  resistance  is  an  important  consideration  for  a  protective  coating. 
However,  the  computational  results  indicate  that  it  would  be  highly  desirable  to  perform 
firing  tests  for  silicon/fiber  coated  fins. 
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Figure  1.  Computational  grid. 


CD 

CD 

CM 

CD 

• 

• 

• 

• 

• 

O 

O 

o 

O 

O 

8 


Figure  2.  In  Boro  tomperatun  font, ours,  aluminum  oxide  coating,  time  ~  5  milliseconds. 
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Figure  5.  In- Bore  temperature  ronton 
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Figure  0.  Free  flight,  temperature  contours,  aluminum  oxide  coating,  time  —  2.0  second;- 
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Figure  14.  Free  flight,  temperature  contours,  silicon/fil 


21 


Figure  1  !>.  Free  flight,  trniprrntim’  contours,  silicon/fil 
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Figure  1G.  Free  flight,  temperature  contours,  silieon/fiher  coating,  time  —  2.5  secomls. 
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Figure  18.  Temprrnt.urr  at.  the  fi 
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APPENDIX  A:  CODE  IMPLEMENTATION  ON  THE  BRL 

SUPERCOMPUTER 


The  heat  transfer  code  is  currently  installed  and  running  on  the  Cray  2  and  Cray 
X-MP/48  Supercomputers  at  BRL.  The  code  is  maintained  using  a  Unix  Makefile.  This 
recompiles  and  links  only  the  code  sections  which  have  been  changed  since  the  last  time  the 
program  was  complied  and  linked.  This  provides  an  efficient  way  to  accurately  maintain 
the  source  code.  The  Makefile  creates  the  executable  code  which  is  then  run  on  the  Cray. 
A  copy  of  this  Makefile  follows: 


#echo 

#echo 

#echo 

#echo 

#echo  The  Subroutines  Are  Being  Compiled  With  Optimization  Turned  On 
#echo 

#echo  iiiiiniiiiiiimmnmnmmimmmmiiiiimmiiMiiimimmiimMmiiiimiiimmimmmi 
#echo 
# 

# 

# 

#geom.h  phyvar.h  solver. h 

#include/geom.h  include/phvvar.h  include/ solver. h 
#include/geom.h  include/phvvar.h  include/solver. h 

# 

& 


2dfin:  nmangfa.o  nequafa.o  nmathfa.o 

segldr  -o  2dfin  nequafa.o  nmangfa.o  nmathfa.o 

# 

# 


echo 
echo 
echo  ’ 
echo 
echo 
echo 
echo  ’ 
echo 
echo  ’ 


LINKING  PROGRAM’ 


OPTIMIZATION  ON’ 
DONE’ 


# 


27 


nmarxgfa.O:  mangfa.f  include/geom.h  include/phvvar.h 
include/sol  ver.h 

echo  '  COMPILING  MANGFA.F' 

/lib/cpp  -P  mangfa.f>nmangfa.f 
cftTT  nmangfa.f 


# 

# 

# 

nequafa.o:  equafa.f  include/geom.h  include/phvvar.h 
include/ solver  .h 

echo  '  COMPILING  EQUAFA.F' 

/lib/cpp  -P  equafa.f>nequafa.f 
cftTT  nequafa.f 


# 

# 

# 

nmathfa.o:  mathfa.f  include/geom.h  include/phvvar.h 
include /sol  ver.h 

echo  ’  COMPILING  MATHFA.F' 

/lib/cpp  -P  mathfa.f>nmathfa.f 
cftTT  nmathfa.f 

Everything  in  the  Makefile  which  is  preceded  with  a  is  a  comment.  The  following 
line  from  the  Makefile  sets  up  the  dependencies  for  the  executable  code. 

2dfin:  nmangfa.o  nequafa.o  nmathfa.o 

segldr  -o  2dfin  nequafa.o  nmangfa.o  nmathfa.o 

W  hat  this  means  is.  the  executable  code  ”2dfin"  depends  on  the  three  object  code  files 
which  follow'  it.  They  are  linked  using  the  Cray  2  linker  ’’segldr.”  All  through  the  Makefile 
you  see  lines  which  are  ” indented”  from  the  first  column.  These  MUST  be  indented  with 
a  ’’tab”  character,  NOT  by  spaces.  This  is  a  restriction  of  the  Makefile,  not  the  operating 
system. 

The  following  line  sets  up  the  dependencies  of  the  ’’nmangfa”  object  code.  This  line 
means  the  object  code  depends  on  the  Fortran  source  code  file  and  the  four  include  files 
which  are  located  in  a  subdirectory  called  ’’include.” 

nmangfa.o:  mangfa.f  include/geom.h  include/phyvar.h 
include/solver.h 

echo  ’  COMPILING  MANGFA.F' 

/lib/cpp  -P  mangfa.f>nmangfa.f 
cft7T  nmangfa.f 
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Due  to  the  length  of  the  code  and  the  numerous  ’’common”  statements,  the  code  is 
set  up  to  use  "include"  statements.  This  allows  changes  to  be  made  in  the  ”.h"  file  and 
the  Makefile  then  causes  these  changes  to  be  made  throughout  the  entire  code.  However, 
the  include  statement  is  not  part  of  standard  Fortran,  so  the  code  is  first  run  through  the 
C-Preprocessor  which  actually  brings  in  the  common  blocks  stored  in  the  .h’  files.  This 

is  done  with  the  ”/lib/cpp  -P . ”  line  of  the  Makefile.  The  output  of  this  preprocessor  is 

redirected  to  the  file  ’’nxxxxx.  f"  which  is  the  Fortran  source  code  with  the  common  blocks 
replacing  the  include  statements. 

To  use  the  include  statements,  the  following  line 

^include”  include /xxxx.h” 


must  be  in  the  Fortran  code,  where  ’’xxxx.h”  is  the  name  of  the  file  to  be  included  and  is 
located  in  a  subdirectory  called  "include."  This  line  goes  in  the  Fortran  code  in  the  place 
you  want  the  ”.h"  file  to  be  inserted.  It  must  start  in  the  FIRST  column,  not  the  seventh 
like  the  Fortran  code  does. 

Due  to  the  length  of  time  necessary  for  some  runs  of  the  code  to  complete,  it  is 
submitted  to  the  batch  queues  of  the  Cray.  This  also  protects  the  run  from  an  unexpected 
system  crash  since  the  code's  point  of  execution  is  quickly  saved  if  the  system  crashes, 
and  its  execution  continues  at  that  point  when  the  system  is  brought  back  up.  Running 
the  code  in  this  way  .uso  permits  you  to  use  the  terminal  for  other  work  while  the  code  is 
running.  To  sr.br  :  the  code  for  execution,  a  short  JCL  file  is  created  which  contains  the 
path  to  the  executable  code  and  any  input  data  files.  The  JCL  file  also  contains  the  name 
of  the  queue  in  which  to  run  the  code,  and  the  maximum  amount  of  CPU  time  to  allocate 
to  it.  There  is  also  a  line  in  the  JCL  file  to  start  and  stop  timing  statistics.  This  makes  it 
easier  to  determine  if  the  code  can  run  in  a  faster  queue  the  next  time  it  is  submitted  since 
thr  queues  are  based  on  the  amount  of  CPU  time  required  and  the  amount  of  memory 
required.  An  example  of  the  JCL  file  follows: 


#  QSUB-q  crayque 

#  QSUB-eo 

#  QSUB-lm  lmw 

#  QSUB-1M  lmw 

#  QSUB-lt  3000 

#  QSUB-1T  3000 

#  QSUB 


#  submit  job  to  the  pipe  queue 

#  standard  error  goes  to  standard  out 

#  establish  per-process  memory  size  limit 

#  establish  per-request  memory  size  limit 

#  establish  per-process  cpu  time  limit  (seconds) 

#  establish  per-request  cpu  time  limit  (seconds) 

#  end  of  QSUB  parameters 


cd  /lfd/sb/eferry/Q3F 
2dfin  <  inp 

The  first  line  of  the  JCL  file  sends  the  job  to  the  queue  ’’crayque.”  This  is  a  ’’pipe” 
queue.  There  are  a  series  of  six  queues  that  the  job  may  run  in.  The  system  looks  at  the 
requirements  of  the  job  and  compares  this  to  the  limitations  of  the  six  queues.  The  job 
is  then  submitted  into  the  most  appropriate  queue.  The  next  line  combines  the  standard 
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error  file  and  the  standard  output  file  into  one,  the  standard  output  file.  The  next  two  lines 
set  up  the  memory  requirements  for  the  code.  Here,  the  requirements  are  1  mega-word. 
The  fifth  and  sixth  lines  in  this  JCL  file  set  the  maximum  time  limit  of  execution  to  3000 
seconds  CPU  time.  The  next  line  supplies  the  path  to  the  code  to  be  run.  The  "2dfin''  is 
the  actual  executable  code,  and  the  ”<  inp’'  redirects  Unix  standard  input  from  the  file 
”inp"  rather  than  from  the  keyboard  of  the  terminal. 

As  the  code  is  executing,  the  real  time  of  flight  is  being  calculated  at  each  time  step. 
When  this  time  just  exceeds  predetermined  values,  a  data  file  (for040.dat)  is  written  out 
with  the  fin’s  grid  geometry  and  the  temperature  at  all  the  grid  points.  This  file  can  be 
plotted  to  the  screen  from  either  of  the  Crays  using  a  program  called  ’’disgrid.”  This  file  can 
also  be  shipped  to  an  Iris  Silicon  Graphics  work  station  using  ftp  (file  transfer  protocol). 
On  the  Iris,  this  file  can  be  graphically  displayed  using  an  Iris  version  of  ’’disgrid.”  This 
’’for040.dat”  file  can  also  be  plotted  using  a  program  called  ”gen”  which  gives  a  smooth 
blending  of  colors  representing  the  temperatures  across  the  grid.  The  "for040.dat’’  files  can 
also  bv.  run  through  another  program  called  "plotc”  which  creates  a  ’’plotc.dat"  file  con¬ 
taining  specific  temperature  contour  levels  specified  within  the  program.  This  ’’plotc.dat" 
file  can  be  plotted  on  the  Crays  using  a  program  called  ’’discont,”  and,  like  the  ’’for040.dat” 
files,  the  ’  plotc.dat"  file  can  be  ftp’d  to  the  Iris  and  be  plotted  by  an  Iris  version  of  ’’dis¬ 
cont.”  The  code  on  the  Iris  which  creates  these  displays  is  written  in  Fortran  and  calls  on 
the  Iris’  extensive  graphics  library  to  visually  display  the  current  state  of  the  fin.  This  can 
be  done  in  one  of  two  ways.  One  way  is  to  view  single  data  files  at  a  time.  This  is  helpful  to 
make  sure  the  code  appears  to  be  running  correctly  before  you  let  it  go  to  its  completion. 
It  is  also  useful  for  creating  hard  copies  of  the  fin  on  a  Tektronix  4692  color  inkjet  printer. 
A  second  method  by  which  the  data  can  be  viewed  is  accomplished  by  reading  in  all  the 
data  files  at  one  time  and  display  them  in  rapid  succession.  This  allows  the  viewer  to  get 
a  feel  as  to  how  the  heat  is  transferring  through  the  fin. 

The  second  method  of  display  uses  a  mode  called  ’’double  buffered  mode”  to  make  the 
successive  images  appear  on  the  screen  more  fluidly.  In  this  mode,  as  one  image  is  being 
viewed,  the  next  image  is  being  drawn  in  a  buffer  ’’behind”  it.  This  image  is  invisible  until 
the  buffers  are  swapped.  When  the  buffers  are  swapped,  the  old  image  gets  pushed  to  the 
back.  This  method  is  very  useful  to  more  easily  see  how  the  fin  changes  from  one  time  to 
the  next.  This  is  especially  useful  to  observe  the  quasi-3d  melting  process. 

The  ’  gen”  program  on  the  Iris  workstation  was  written  by  Harry  Dwyer  using  GSS  on 
an  IBM  compatible  PC,  and  was  ported  to  the  Iris  by  Earl  N.  Ferry,  Jr.  The  program  uses 
a  method  of  shading  on  the  Iris  known  as  Gouraud  shading  to  achieve  smooth  contours. 
This  is  accomplished  by  creating  a  color  map  and  interpolating  the  temperatures  of  the 
grid  onto  the  indices  of  the  color  map.  These  indices  are  then  used  as  color  shades,  for 
the  fill  routine.  This  method  of  coloring  the  grid  creates  a  smooth  blend  of  colors  for  the 
different  temperatures. 

The  heat  transfer  code  also  creates  a  file  called  ’’tdat.dat.”  This  file  contains  the 
temperatures  of  specific  grid  points  at  each  time  interval.  This  file  can  be  plotted  on  the 
Crays  using  a  program  called  ’’dta.x.”  This  program  plots  this  file  as  temperature  vs.  time 
for  the  selected  points. 
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The  programs  which  are  ran  on  the  Cravs  to  plot  the  grids  and  contours  are  written 
in  Fortran  and  currently  use  DISSPLA.  These  programs  also  present  another  hard  copy 
option.  The  "disgrid"  and  ”discontv  programs  have  the  option  to  create  a  fii<=>  which  can 
be  ftp’d  to  the  Yax  and  be  printed  on  a  QMS  Laser  printer. 
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Copies  Organization 

1  Office  of  the  Secretary  of  Defense 
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Director,  Live  Fire  Testing 
ATTN:  James  F.  0  Bryon 
Washington,  DC  20301-3110 


No  of 

Copies  Organization 
1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
Ames  Research  Center 
Moffett  Field,  CA  94035-1099 


2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
.Alexandria,  VA  223CM-6145 

1  HQDA  (SARD-TR) 

WASH  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 


1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

1  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Library) 
Warren,  MI  48397-5000 

1  Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002-5502 


1  Commander 

US  Armv  Laboratory  Command 
ATTN:  AMSLC-DL 
Adclphi,  MD  20783-1145 


(Class,  only)  ]  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 
Fort  Benning,  GA  31905-5660 


Commander 
US  Army,  ARDEC 
ATTN:  SMCAR-1MI-I 
Picatinny  Arsenal,  NJ  07806-5000 


(tnciM*.  only)  ]  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 


2  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR  TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-DACL 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 


1  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

1  Cdr,  USATECOM 

ATTN:  AMSTE-TD 

3  Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

1  Dir,  VLAMO 

ATTN:  AMSLC-VL-D 
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No  of 

Copies  Organization 


4  Commander 

Armament  RD&E  Center 
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Chung 

Kahn 
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3  Commander 

Armament  RD&E  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-CCL-CA, 

Hirlinger 

O’Niell 

Miller 
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1  Commander 
Armament  RD&E  Center 
US  Army  AMCCOM 

ATTN:  SMCAR-FSP-A(l)/Bischer 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  DK20  (Clare) 

Code  DK20  (Moore) 
Dahlgren,  VA  22448-5000 

2  Director 

Sandia  National  Laboratories 
ATTN:  Dr.  W.  Oberkampf 
Dr.  F.  Blottner 
Division  1636 
P.O.  Box  5800 
Albuquerque,  NM  87185 

1  Massachusetts  Institute  of  Technology 
ATTN:  Technical  Library 

77  Massachusetts  Avenue 
Cambridge,  MA  02139 

2  Commander 

David  W.  Taylor  Naval  Ship 
ATTN:  Dr.  S.  de  los  Santos 
Mr.  Stanley  Gottlieb 
Bethesda,  MD  20084-5000 

1  Commandant 
IISAFAS 

ATTN:  ATSF-TSM-CN 
Fort  Sill,  OK  73503-5600 

1  University  of  California,  Davis 

Department  of  Mechanical  Engineering 
College  of  Engineering 
ATTN:  Prof.  Harry  A.  Dwyer 
Davis,  CA  95616 


2  Ford  Aerospace  and  Communications 
Corporation 
Aeronautronic  Division 
ATTN:  Charles  White 
Budd  Blair 
Ford  Road 

Newpoint  Beach,  CA  92658 

1  United  States  Military  Academy 

Department  of  Civil  &  Mechanical  Engineering 
ATTN:  LTC  Andrew  L.  Dull 
West  Point,  NY  109% 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1  not  vt  W  BRL-MR-3831  APR1L  1990 

1 .  BRL  Report  Number _ _ Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  cic0  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to'  organization,  technical  content,  format,  etc.) _ 


Name 


CURRENT  Organization 
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Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 
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Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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